home *** CD-ROM | disk | FTP | other *** search
/ Just Call Me Internet / Just Call Me Internet.iso / com / computer / casio_st / prog_fx / oisac2 / l_s6420.txt < prev    next >
Text File  |  1994-03-01  |  575b  |  15 lines
  1. 10 GOSUB400:IFp>1 ORp=<0 THEN*
  2. 15 IFp=1 THENx=0:GOTO80
  3. 20 IFn=1 THENq=p:p=q/2:GOSUB"LIB0:S6410":GOSUB400:x=y*y:p=q:GOTO80
  4. 30 IFn=2 THENx=-2*LNp:GOTO80
  5. 40 GOSUB"LIB0:S6410":GOSUB400:x=y:y=x*x:x=n+SQR(2*n)*x+2/3*(y-1)+x*(y-7)/9/SQR(2*n)-2/405/n*(y*(3*y+7)-16):IFx<0 THENx=0
  6. 50 IFn>100 THEN80
  7. 60 b=x:GOSUB100:IFn<31 THENIFABS(b-x)>1E-4 THEN60
  8. 80 y=x:RETURN
  9. 100 q=p:GOSUB"LIB0:S6220":GOSUB400:p=y
  10. 110 z=EXP(((n-1)*LN(x/n)-LN(4*PI*x)+n-x-1/n/6)/2)
  11. 120 x=x+(p-q)/z:p=q:x=ROUND(x,-6):RETURN
  12. 400 ONERRORGOTO500:RETURN
  13. 500 IFERR=1 THENCLS:ONERRORGOTO0
  14. 510 sr=1:RESUME80
  15.